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ABSTRACT 


This  technical  note  presents  a  Kalman  filtering  approach  that  is  used  for  filtering  narrow- 
band  interferers  out  of  direct  sequence  spread  spectrum  signals.  The  approach  is  based 
on  the  digital  phase-locked  loop  Kalman  filter  and  is  close  to  optimum  in  so  far  as  de¬ 
modulating  an  FM-type  of  interferer.  Because  the  interference  is  assumed  to  be  much 
stronger  than  either  the  signal  or  noise,  the  Kalman  filter  locks  onto  the  interference 
and  produces  estimates  of  its  phase  and  envelope.  The  algorithm  is  compared,  through 
computer  simulation,  with  the  recursive  least  squares  lattice  algorithm  for  the  case  of  a 
swept  tone  interferer.  Examples  of  the  phase-  and  envelope-tracking  capabilities  of  the 
algorithm  are  presented,  followed  by  bit  error  rate  curves  for  the  case  of  an  interferer  of 
bandwidth  equal  to  20%  of  the  chip  rate  of  the  spread  spectrum  signal.  The  results  show 
that  with  proper  selection  of  the  Kalman  filter  parameters,  the  Kalman  filter  excisor  can 
outperform  the  recursive  least  squares  algorithm  by  2  dB.  It  is  suggested  that  performance 
could  perhaps  be  improved  further  by  using  a  higher  order  Kalman  filter  to  remove  the 
residual  phase-tracking  error  produced  by  the  second  order  Kalman  filter. 

RESUME 

*  Ce  rapport  contient  une  description  d’un  filtrage  de  Kalman  utilise  pour  enlever  des  in¬ 

terferences  a  bande  passante  etroite  d’un  signal  a  spectre  etale  par  une  sequence  directe. 
La  methode  repose  sur  un  filtre  de  Kalman  a  verrouillage  numerique  de  phase  et  est 
presque  optimale  pour  la  demodulation  d ’interferences  de  type  FM.  Le  filtre  de  Kalman 
se  verrouille  sur  I’interference  et  produit  une  estimation  de  sa  phase  et  de  son  envelope 
parce  que  I’on  suppose  que  I’interference  et  le  bruit  sent  beaucoup  plus  forts  que  le  sig¬ 
nal.  Une  simulation  par  ordinateur  est  utilise  pour  comparer  ce  filtre  a  la  methode  des 
treillis  des  moindres  carres  recursifs  pour  une  interference  a  balayage  de  frequence.  On 
presente  des  exemples  de  poursuites  de  phase  et  d’enveloppe  pax  I’algorithme  ainsi  que 
des  courbes  de  taux  d’erreur  des  bits  pour  une  interference  dont  la  bande  passante  est 
de  20%  du  debit  numerique  du  signal  a  spectre  etale.  II  est  demontre  qu’une  selection 
judicieuse  des  parametres  du  filtre  de  Kalman  permet  une  excision  meilleure  que  celle 
produite  par  I’algorithme  recursif  des  moindres  carres  par  2  dB.  II  est  suggere  qu’une 
amelioration  supplementaire  pourrait  etre  obtenue  par  I’utilisation  d’un  filtre  de  Kalman 
d’ordre  superieur  qui  enleverait  ce  qui  reste  de  I’erreur  de  phase  d’un  filtre  de  Kalman  de 

,  second  ordre. 
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EXECUTIVE  SUMMARY 


This  technical  note  presents  a  Kalman  filtering  approach  that  is  used  for  filtering  narrow- 
band  interferers  out  of  direct  sequence  spread  spectrum  signals.  These  signals  are  used 
extensively  in  military  communication  systems.  The  technique  described  herein  applies 
equally  to  both  Electronic  Support  Measures  (ESM)  systems  and  direct  sequence  spread 
spectrum  communication  sytems.  In  the  former  application,  the  ESM  system  may  be 
attempting  to  intercept  the  spread  spectrum  signal,  but  the  narrowband  interference  may 
be  hampering  this  effort.  In  the  latter  application,  the  spread  spectrum  communication 
system  may  require  additional  assistance  to  suppress  the  interference.  Since  the  open 
literature  has  been  devoted  to  this  latter  case,  the  material  presented  here  focuses  on  this 
application. 

One  of  the  attributes  of  direct  sequence  spread  spectrum  communication  systems 
is  their  ability  to  combat  interference  or  intentional  jamming  by  virtue  of  the  system’s 
processing  gain  inherent  in  the  spreading  and  despreading  process.  The  interference  can  be 
attenuated  by  a  factor  up  to  this  processing  gain.  In  some  cases  the  gain  is  insufficient  to 
effectively  suppress  the  interferer,  leading  to  a  significant  degradation  in  communications 
manifested  by  a  sudden  increase  in  bit  error  rate.  If  the  ratio  of  interference  bandwidth 
to  spread  spectrum  bandwidth  is  small,  the  interference  can  be  filtered  out  to  enhance 
system  performance.  However,  this  is  at  the  expense  of  introducing  some  distortion  onto 
the  signal.  This  process  of  filtering  is  sometimes  referred  to  as  interference  excision. 

The  Kalman  filtering  approach  is  based  on  the  digital  phase-locked  loop  Kalman 
filter  and  is  close  to  optimum  in  so  far  as  demodulating  an  FM-type  of  interferer.  Because 
the  interference  is  <issumed  to  be  much  stronger  than  either  the  signal  or  noise,  the  Kalman 
filter  locks  onto  the  interference  and  produces  an  estimate  of  the  phase  and  envelope  of 
the  interference. 

The  algorithm  is  compared,  through  computer  simulation,  with  the  recursive 
least  squares  lattice  algorithm  for  the  case  of  a  swept  tone  interferer.  E.xamples  of  the 
phase-  and  envelope-tracking  capabilties  of  the  algorithm  are  presented,  followed  by  bit 
error  rate  curves  for  the  case  of  an  interferer  of  bandwidth  equal  to  20%  of  the  chip  rate 
of  the  spread  spectrum  signal.  The  results  show  that  with  proper  selection  of  the  Kalman 
filter  parameters,  the  Kalman  filter  excisor  can  outperform  the  recursive  least  squares 
algorithm  by  2  dB  and  that  performance  could  perhaps  be  improved  further  by  using  a 
higher  order  Kalman  filter  to  remove  the  residual  phase-tracking  error  produced  by  the 
second  order  Kalman  filter  used  in  this  study. 
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1.0  INTRODUCTION 


Direct  sequence  spread  spectrum  communication  systems  have  an  inherent  processing 
gain  which  can  reduce  the  effects  of  jammers  or  unintentional  interference.  When  these 
intruding  signals  have  a  power  advantage  over  the  spread  spectrum  system,  a  severe  degra¬ 
dation  in  communications  results.  However,  communications  can  be  enhanced  somewhat 
by  filtering  the  interference,  particularly  if  its  bandwidth  is  significantly  less  than  the 
bandwidth  of  the  spread  spectrum  signal. 

A  host  of  digital  filtering  algorithms  based  on  time-modelling  concepts  have 
been  developed  over  the  years  to  address  this  problem  [1,  2,  3].  The  received  signal, 
noise  and  interference  are  applied,  for  example,  to  a  filter  matched  to  a  chip  and  the 
output  is  sampled  at  the  chip  rate.  If  one  assumes  that  the  resulting  sampled  data  can 
be  modelled  as  an  autoregressive  (AR)  process,  then  the  AR  coefficients  or  their  related 
lattice  reflection  coefficients  can  be  determined  from,  say,  the  maximum  entropy  algorithm 
[4],  stochastic  gradient  algorithms  [5],  and  block  and  recursive  least  squares  algorithms 
[6,  7].  Estimates  of  these  coefficients  by  the  above  algorithms  lead  directly  to  transversal 
or  lattice  filter  structures,  both  of  which  act  as  whitening  filters  for  the  assumed  AR 
process.  This  approach  to  suppressing  interference  (or  excision  as  it  is  sometimes  referred 
to)  is  possible  because  of  the  non-coherency  of  the  signal  and  noise  samples,  and  the  more 
coherent  interference  samples  resulting  from  the  latter’s  narrowband  nature. 

A  different  time-modelling  approach  to  interference  excision  was  presented  in 
[8].  There,  a  function  related  to  the  interference  is  assumed  to  have  been  generated  by 
a  second  order  state-space  system  to  which  the  extended  Kalman  filter  equations  are 
applied.  This  results  in  what  has  been  termed  the  digital  phase-locked  loop  (DPLL) 
[9,  10].  This  approach  to  interference  suppression  is  close  to  optimum  for  the  suppression 
of  constant  envelope  or  FM-types  of  interferers  and,  therefore,  can  potentially  provide 
better  performance  than  the  recursive  least  squares  (RLS)  algorithm  which  was  analyzed 
for  the  swept  tone  in  [11].  The  objective  of  this  technical  note  is  to  present  a  performance 
comparison  between  the  DPLL  and  RLS  algorithms  for  the  swept  tone  interferer. 

The  outline  of  this  technical  note  is  as  follows.  Section  2.0  describes  the  spread 
spectrum  communication  model  and  interference  estimator  based  on  the  Kalman  filter. 
Section  3.0  summarizes  the  interference  estimator.  Section  4.0  presents  simulation  results 
for  the  case  of  a  swept  tone  interferer  corrupting  the  spread  spectrum  signal.  The  Kalman 
filtering  approach  is  compared  to  the  lattice  RLS  algorithm  previously  discussed  [llj. 
Finally,  Section  5.0  concludes  the  technical  note,  suggesting  areas  for  further  research. 


2.0  COMMUNICATIONS  MODEL 


The  basic  elements  of  the  BPSK  PN  spread  spectrum  receiving  system  are  shown  in 
Fig.  1.  The  received  waveform  r(t),  consisting  of  a  spread  spectrum  signal,  additive 


PN 

Correlator 


Figure  1:  Spread  spectrum  communications  model. 

white  Gaussian  noise,  and  narrowband  interference  is  applied  to  a  bandpass  filter  with 
the  transfer  function  Hbp(f),  whose  output  is  defined  as 

u{t)  =  s{t)  +  n{t)  +  i{t).  (1) 

The  bandpass  filter  for  the  application  considered  here,  is  assumed  to  be  a  filter 

matched  to  a  chip  and  centered  at  the  carrier  angular  frequency  ujq  of  the  spread  spectrum 
signal.  The  spread  spectrum  signal  is  defined  as 

s(t)  =  a(t)cos(a;oO  (2) 

where 

a{t)  =  J2Dkbkit-kTb).  (3) 

k 

In  Eq.  (3),  Dk  is  a  sequence  of  data  bits  of  amplitude  (±1)  and  duration  Tj,  seconds,  and 
bk{t  —  kTb)  is  the  PN  sequence  pattern  for  the  bit,  i.e., 

bk{t)  =  J2^kjqit-jTc)  (4) 

i=i 
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with  L  being  the  number  of  pseudo  random  chips  per  bit,  or  the  processing  gain,  is 
the  code  sequence  for  the  bit,  and  q{t)  represents  the  basic  chip  pulse  of  energy  Ec- 

The  noise  n{t)  is  Gaussian  and  hcis  a  power  spectral  density 

SM)  =  yI%(/)P,  (5) 

where  No/2  is  the  power  spectral  density  of  the  assumed  white  Gaussian  noise  from  the 
channel.  The  band  of  interference  is  defined  as 

i{t)  =  I{t)  cos{ujot  +  <?(<))  (6) 

where  /(<)  is  the  interference  amplitude  and  6{i)  is  the  phase  modulation.  It  has  been 
assumed  that  the  effect  of  the  bandpass  filter  Hbp{f)  is  negligible  on  the  interference  i{t). 
This  would  not  be  true  in  reality,  since  the  bandpass  filter  will  induce  some  amplitude 
modulation  on  the  constant  envelope  interference.  The  rate  of  change  of  the  envelope 
would  therefore  depend  on  the  rate  of  change  of  the  instantaneous  frequency,  and  the 
amplitude  deviation  would  depend  on  the  frequency  offset  from  uq  and  the  interference 
bandwidth.  The  fact  that  this  report  considers  the  amplitude  to  be  of  constant  envelope 
at  the  output  of  the  bandpass  filter  corresponds  to  the  best  case  situation. 

Referring  to  Fig.  1,  the  output  u{t)  of  the  bandpass  filter  Hbp(f)  is  bandpass 
sampled  and  applied  to  a  limiter/bandpass  filter  and  interference  estimator. 

Consider  the  bandpass  sampler  first.  The  analog  signal  u{t)  from  Eq.  (1)  is 
sampled  at  /,  =  2Rc  {pifs  =  a;o/2;r  +  /2c/2  for  some  integer  m)  where  Rc  is  the  chip  rate 
of  the  spread  spectrum  signal.  The  resultant  sampled  signal  is,  therefore, 

Un  =  -Sn  +  n„  +  l„,  (7) 

where  s„  consists  of  the  sequence  {. . .  ,0, (— l)"a„,0,  (  — •  •  •}  where  the  a„  are 
of  energy  Ec  and  coded  according  to  CkjDk  for  the  chip  in  the  transmitted  bit. 
Tin  ^  are  uncorrelated  Gaussian  noise  samples  of  variance  <7^,^  =  Ec{No/2),  and  is  the 
sampled  version  of  Eq.  (6).  The  samples  Un  are  applied  to  the  interference  estimator  and 
interference  canceller. 

The  interference  estimator  produces  an  estimate  of  the  interference  which  is 
removed  from  the  sampled  input  Un-  The  output  of  the  summer  is  decimated  and  sign- 
reversed,  resulting  in  a  baseband  error  signal  e„.  This  error  signal  is  correlated  with  the 

'Coherent  bandpass  sampling  has  been  assumed,  so  that  the  in-phase  component  is 
(«i,n/V2)  co8(n»'/2)  and  the  quadrature  component  is  (n2,„/\/2)  sin(njr/2) 
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PN  sequence.  The  output  of  the  correlator  is  integrated  and  bit-detected. 


3.0  INTERFERENCE  ESTIMATOR 


Consider  now  the  branch  containing  the  limiter.  Here,  u„  is  applied  to  a  limiter/bandpass 
filter.  The  input  to  the  limiter  referenced  to  the  interference  is  redefined  as 


«n  =  \J[h  +  n[  „  -I-  -1-  [n'2,„  +  cosfcjon  +  0^  +  <i>u.n] 


(S) 


where 


flirctstri 


«2.n/>/2  +  Oj.n 


(9) 


In  +  ni,„/\/2  +  ai,„ 

is  a  noise- like  phase  fluctuation  on  the  interferer’s  phase  and  is  due  to  the  noise  and 
spread  spectrum  signal.  The  terms  nj  and  are  in-phase  and  quadrature 

components  of  the  noise  and  spread  spectrum  signal  with  respect  to  the  interference  phase 

dn. 

The  output  of  the  limiter/bandpass  filter  is  [12] 


iA' 

i>Pn  =  - cosfwon  +  +  <t>u,n]^ 

TT 


(10) 


where  A'  is  the  limiter’s  output  level.  This  signal  is  redefined  as  (letting  A'  =  ■\/2x/4  for 
convenience) 

hpn  =  \/2  cosjwon  A  9^  A  <i>u,n] •  (11) 

It  should  be  noted  that  for  large  interference-to-noise  ratios  in  which  the  interference  is 
of  constant  envelope,  in  Eq.  (11)  is  approximately  Gaussian  [12].  The  sampled  signal 
in  Eq.  (11)  is  what  is  processed  by  the  Kalman  filter,  which  estimates  the  phase  O-n  of  the 
interference.  The  phase  estimate  is  denoted  as 

The  interference  estimator  shown  in  Fig.  1  is  detailed  in  Fig.  2.  The  Kalman 
algorithm  [8]  produces  a  signal 


yn  =  V2  cos(u;on  A  K]n-\)  (12) 

which  can  be  used  as  a  basis  for  estimating  a  sampled  version  of  the  envelope  of  i(t)  (i.e., 
In  in  Eq.  (6)).  The  output  of  the  first  multiplier  is  a  baseband  term  and  is,  using  Eq.  (7), 
excluding  the  \/2  factor. 

In  ~  In  COs(fln  ^n|n— l) 


4 


+  n„  cos(u;an  +  +  a„  cos(^„|„_i) 

+  a„cos(2a;on  + 


(i;}) 

Equation  (13)  consists  of  four  terms:  the  first  term  is  related  to  the  desired  envelope  of 


Figure  2:  Block  diagram  of  the  interference  estimator. 

the  interference;  the  second  is  approximately  Gaussian  baseband  noise  [13];  and  the  third 
and  fourth  terms  are  noise-like  terms  emanating  from  the  spread  spectrum  signal.  The 
fourth  term,  because  of  the  sampling  rate  conditions  discussed  in  [9],  is  essentially  filtered 
out  by  the  low  pass  filter  of  bandwidth  Bipjr  <  0.50  Hz  and,  therefore,  will  be  ignored  in 
the  baseband  simulations  to  be  discussed  in  the  next  section.  The  term  /^/\/2  is  filtered, 
resulting  in  the  estimate  of  the  interference  envelope,  /„/\/2.  Combining  this  with  y„  in 
Eq.  (12)  and  shown  in  Fig.  2  yields,  the  estimate  of  the  interference 

in  =  /n  cos(u;on  +  0„|„_i),  (14) 

which  is  subtracted  from  «„  as  illustrated  in  Fig.  1. 

4.0  SWEPT  TONE  PERFORMANCE 

The  objective  of  this  section  is  to  compare  the  performance  of  the  Kalman  filter  excisor 
to  the  RLS  algorithm  [11]  for  the  case  of  a  swept  tone  interferer.  For  this  comparison, 
the  Kalman  filter  is  not  adaptive;  an  adaptive  architecture  will  be  discussed  in  another 
technical  note.  It  is  intended  here  to  illustrate  the  potential  performance  of  the  Kalman 
filtering  approach.  The  modified  Kalman  algorithm  [8]  is  the  one  that  will  be  used. 

From  [8]  the  performance  of  the  Kalman  filter  is  governed  by  three  parameters. 
These  are  the  frequency  deviation  constant  d,  the  variance  of  the  observation  noise  cr^  = 
Nobsl^  assumed  in  the  model  where  Nobal^  is  the  two-sided  power  spectral  density  of  the 
observation  noise,  and  the  3  dB  bandwidth,  a/,  of  a  low  pass  filter  for  filtering  white 
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Gaussian  noise  to  produce  the  interferer’s  modulating  signal.  It  is  important  to  note 
that  a/,  which  is  related  to  7  =  /^/ay,  the  ratio  of  the  sampling  rate  to  this  bandwidth, 
establishes  a  family  of  performance  curves  such  as  those  in  [9].  Thus,  for  a  given  aj  and 
the  approach  to  be  taken  will  be  to  determine  that  value  of  d  which  minimizes  the 
mean-squared  value  of  the  residual  interference  in  front  of  the  PN  correlator  in  Fig.  1. 
It  is  assumed  here  that  for  small  enough  Bipp  in  Fig.  2,  the  main  source  for  residual 
interference  will  be  the  phase-noise.  Therefore,  using  the  minimum  mean-squared  value 
of  the  residual  interference  as  the  criterion  for  optimality,  the  value  of  d  determined  would 
provide  as  close  to  the  optimum  performance  achievable.  Finally,  the  variance  parameter 
for  the  observation  noise,  crl,  will  be  set  according  to  the  noise  level  in  the  experiment. 
This  noise  level,  as  shown  below,  is  a  function  of  the  power  spectral  densitie.s  of  the 
channel  noise  and  spread  spectrum  signal. 

For  this  experiment,  the  interference  is  modelled  as  a  frequency  modulated  wave¬ 
form  of  the  type 


i{t)  =  I  cos{uQt  +  6{t)) 

=  I  cos{ujQt  + ‘ZirS ft  +  d  f  m{T)dT  +  ip)  (15) 

Jo 

where  6f  is  an  offset  frequency  from  the  spread  spectrum  signal  carrier,  m(r)  is  the 
modulating  signal  of  the  interferer,  d  is  the  frequency  deviation  constant  of  the  modulator 
in  radians/second/ volt  [14],  and  <fi  is  an  initial  random  ph^lse.  The  swept  tone  with  sweep 
rates  of  ±0.001  Hz/sec.  was  modelled  using  Eq.  (15).  A  section  of  the  instantaneous 
frequency  ^  6{t)/2v  from  Eq.  (15)  is  shown  in  Fig.  3,  sweeping  between  0.1  and  0.3  Hz, 
with  Sf  =  0.2  Hz  and  d  =  0.628  radians/second/volt.  With  this  sweep  rate  and  the 
values  of  Bipp  used  in  Fig.  2,  it  is  assumed  that  the  envelope  variation  at  the  output  of 
the  bandpass  filter  is  slow  enough  to  render  the  constant  envelope  assumption  valid.  The 
peak  frequency  deviation  A/  for  this  case  was  0.1  Hz  and  the  interference-to-t'  ■  '>  ;  itio 
was  set  at  20  dB. 

Concerning  the  other  parameters  of  the  experiment,  the  processing  gain  was  set 
at  L  =  20  and  the  energy/chip  w^ls  Ec  =  I  unit.  For  the  Kalman  filtering  approach,  the 
excisor  used  is  the  one  in  Fig.  1,  with  the  interference  estimator  cis  shown  in  Fig.  2.  In  the 
latter  figure,  the  low  pass  filter  is  of  the  order  Butterworth  type  of  bandwidth  B^^pp, 
and  is  a  variable  in  the  simulations.  For  the  RLS  excisor,  one  can  refer  to  [llj. 

^In  the  simulations  that  follow,  the  carrier  frequency  /o  has  been  assumed  to  be  zero,  i.e.,  baseband 
models  have  been  used. 
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The  experiment  consists  of  two  parts.  The  first  part  is  concerned  witli  deter¬ 
mining  from  simulation  the  optimum  value  of  d  for  a  particular  a*.  The  second  part  is 
concerned  with  obtaining  the  bit  error  rate  curve  for  the  swept  tone  interferer  using  tliis 
optimum  value  of  d,  while  modifying  al  in  relation  to  the  particular  Eb! No- 


4.1  Optimum  d  Determination 

For  fixed  aj  and  it  suffices  to  determine  the  value  of  the  power  spectral  density  of  the 
observation  noise  =  Nobal^  to  be  used  in  the  Kalman  filter  before  varying  d  to  obtain 
the  optimum  value.  This  is  accomplished  as  follows. 

The  interference-to-noise  ratio  INR  in  the  equivalent  noise  bandwidth  of  the 
filter  of  3  dB  bandwidth  a/  in  the  Kalman  filter  interference  model  in  [9]  is 


INR  = 


2 

TCCtj  Nobs 


(16) 


From  Eq.  (9),  the  phase-noise  at  the  output  of  the  limiter  and  bandpass  filter  was  defined 
as 


<f>u,n  =  arctan 


(17) 


where  /  has  replaced  since  a  constant  envelope  interferer  is  being  considered.  In 
Eq.  (17),  using  the  definition  for  the  signal  and  noise  from  Eq.  (7),  the  noise  and  spread 
spectrum  signal  referenced  to  the  phase  of  the  interference  become 


=  "l.nCOs(0„)  -f  n2,nSin(0„) 


”2,n  =  ”2,n  cos(^„)  -  ni.„  sin(^„) 
«l.n  =  anCOs(0n) 

02  „  =  -ansin(^„). 

For  I  {sn  +  rin),  Eq.  (17)  simplifies  to 


4>u,n  « 


The  approximate  mean-squared  value  of  the  phase-noise  is,  therefore, 


^{<n}  «  (iVo/2  +  0.5)//2 


(18) 

(19) 

(20) 
(21) 


(22) 
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=  Nobs/2 


(23) 


where  the  0.5  in  Eq.  (23)  refers  to  the  mean-squared  value  of  the  signal  term  in  Efj.  (22). 
It  can  be  seen  from  Eq.  (23)  that  the  variance  of  the  observation  noise,  cr^,  in  the  Kalman 
filter  model  is  a  function  of  the  channel  noise  power  spectral  density  No,  the  power  in  the 
spread  spectrum  signal,  and  the  power  in  the  interference. 

The  kind  of  performance  that  can  be  achieved  with  the  Kalman  filter  using 
particular  values  for  <tI  and  aj  while  varying  d  is  considered  next.  For  this  example  No 
(and  therefore  al  from  Eq.  (23))  was  determined  from  Eb/No  =  12  dB,  q/  was  set  to  a 
value  of  0.0003125  Hz,  which  corresponds  to  the  value  used  in  [9],  d  ranged  from  0.1  to 
2.0  rad./sec./volt,  and  the  bandwidth  of  the  low  pass  filter  designated  as  Blpf  in  Fig.  2 
was  either  0.10  or  0.025  Hz. 

Figure  4  shows  estimates  of  the  degree  of  interference  suppression  at  the  output 
of  the  Kalman  filter  excisor,  where  suppression  is  defined  as 


5  =  101og(AiyP.) 


(24) 


where  Ain  =  in  —  in  is  the  residual  interference,  and  Pi  is  the  power  in  the  interference. 
The  number  of  samples  used  in  determining  the  estimate  of  <5  in  Eq.  (24)  was  20,000. 
Observe  that  the  “optimum”  value  for  d  is  %  0.20  rad./sec./volt. 

From  Eq.  (24)  and  the  results  in  Fig.  4,  the  efficiency  of  the  excisor,  //,  defined 
in  Eq.  (Ola)  in  [11]  can  be  determined  from 


1 

PO  +  S 


(25) 


where,  for  the  caise  considered  here,  <7*,^  =  Nq/2  and  po  =  Pi.  This  can  be  seen  by  noting 
that,  for  the  linear  prediction  approach  for  which  this  equation  was  developed  [11],  the 
overall  noise  term  in  Eq.  (Ola)  in  [11]  can  also  be  expressed  as 


^2  _ 

-  <^nse 


(26) 


It  should  be  noted,  however,  that  Eq.  (25)  for  the  Kalman  filter  excisor  is  only  approx¬ 
imate,  with  the  accuracy  improving  as  the  bandwidth  of  the  lowpass  filter  in  Fig.  2 
becomes  significantly  less  than  the  sampling  rate  /,.  Under  these  latter  conditions,  the 
cross-correlation  between  the  residual  interference  and  the  signal  and  noise  components 
(which  is  zero  for  the  linear  prediction  case)  at  the  output  of  the  excisor  becomes  small. 

As  an  example  of  the  amount  of  efficiency  that  can  be  achieved  using  Eq.  (26), 
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for  EfifNo  =  12  dB,  I/S  =  20  dB  and  S  =  —22  dB  at  d  =  0.20  rad. /sec. /volt  in  Fig.  1, 
7/  =  19  dB. 


S,  dB 


Figure  4:  Suppression  level  S  from  Eq.  (24)  for  several  values  of  d,  and  B^pf  =  0.1  Hz. 


An  example  of  the  phase-tracking  capability  of  the  Kalman  filter  for  the  swept 
tone  interferer  is  illustrated  in  Fig.  5,  using  the  optimum  value  of  d  =  0.20  rad. /sec. /volt 
from  Fig.  4.  There  is  a  slight  offset  of  the  pheise  estimate  0„|„_i  from  the  exact  phcise 
On-  This  offset  is  27r  radians  and  is  due  to  initial  cycle  slippage  of  the  Kalman  filter. 
Furthermore,  because  of  the  nonlinearly  varying  form  of  the  phase  (it  has  a  “sinusoidal” 
feature  to  it  superimposed  on  a  ramp  of  slope  27r(0.2)T’,  radians/sec.  —  refer  to  Eq.  (15) 
and  the  modulating  signal  in  Fig.  3),  the  estimate  moves  closer  and  further  away  from  the 
exact  phase.  This  phenomenon  is  better  illustrated  in  Fig.  6  which  shows  the  phase  error 
—  ^n|n-i>  with  the  initial  27r  offset  removed,  and  is  due  to  the  inability  of  the  modified 
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Kalman  filter  to  track  nonlinear  phase  functions  without  some  residual  phase  error. 

The  envelope-tracking  capability  of  the  interference  estimator  of  Fig.  2  is  illus¬ 
trated  in  Fig.  7  for  two  cases:  Blpf  =  0.10  and  0.025  Hz.  Observe  that  by  using  a 
smaller  value  of  B^pp,  the  envelope  noise  has  decrecised,  which  is  to  be  expected,  but  at 
the  expense  of  a  longer  time  constant. 

Finally,  the  power  spectra  for  the  interference  and  its  estimate  are  illustrated 
in  Fig.  8.  The  accuracy  of  the  estimate  is  excellent  in  the  region  of  the  spectrum  where 
most  of  the  interference  power  is  located;  beyond  this  region  one  can  see  the  inaccuracies 
imposed  by  the  noise  due  to  both  envelope  and  phase-noise. 


Figure  5:  An  example  of  the  phase-tracking  ability  of  the  Kalman  filter  for  d  =  0.20 
rad./sec./volt,  auid  Bipf  =  0.10  Hz. 
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Figure  7:  An  example  of  the  envelope  estimator  for  d  =  0.20  rad. /sec. /volt,  and  Bipp  = 
0.10  and  0.025  Hz, 
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Frequency  (Hz) 


Figure  8:  Power  spectrum  of  the  interference  and  its  estimate  for  d  =  0.2  rad. /sec. /volt, 
and  BhPF  =  0.1  Hz. 


4.2  Bit  Error  Rate  Performance 


In  this  section  the  bit  error  rate  performance  of  the  Kalman  filter  is  compared  with  the 
RLS  algorithm  discussed  in  [11].  For  the  simulation,  a  total  of  50,000  bits  were  transmiltc'd 
(1,000,000  chips),  with  EblNo  ranging  fromO  to  12  dB.  The  Kalman  fdter  parameters  were 
set  up  as  follows:  ctj  =  0.0003125  Hz,  d  —  0.20  rad  ./sec. /volt  (i.e.,  the  optimum  value 
from  Fig.  4),  and  a  range  of  values  for  <t^,  which  were  calculated  from  Eq.  (23)  for  the 
different  values  of  EblNo.  Finally,  the  lowpass  filter  Bipf  in  Fig.  2  used  to  estimate  the 
envelope  /  of  the  interference,  was  0.10,  0.05  and  0.025  Hz.  The  parameters  of  the  RLS 
algorithm  were  the  same  as  those  used  in  [11]. 

The  results  are  illustrated  in  Fig.  9.  They  indicate  that,  for  the  values  of  Bcpp 
used,  the  Kalman  filter  performs  better  than  the  RLS  algorithm  with  its  matched  filter. 
In  fact,  the  improvement  in  performance  gets  progressively  better  as  Blpf  decreases. 
However,  there  is  little  noticeable  difference  between  Blpf  =  0.05  Hz  and  Blpf  =  0.025 
Hz,  suggesting  that  the  residual  pha.se  error  and  phase-noise,  both  from  the  Kalman 
filter,  as  opposed  to  the  additive  noise,  and,  to  some  degree,  signal  distortion,  are  the 
main  contributors  to  the  lower  bound  in  performance.  Considering  the  nonlinear  phase 
function  which  the  Kalman  filter  is  attempting  track,  it  would  be  interesting  to  see  if  a 
higher  order  Kalman  filter  could  provide  even  better  bit  error  rate  performance. 

5.0  CONCLUSIONS 

This  technical  note  presented  some  of  the  performance  capabilities  of  the  Kalman  filter 
excisor  as  applied  to  a  swept  tone  interferer  with  interference-to-signal  ratio  equal  to 
20  dB.  The  Kalman  filter  technique  is  close  to  optimum  for  this  type  of  interference  as 
compared  to  the  RLS  excisor,  although  the  latter  technique  can  produce  good  results  if 
the  whitening  filter  is  followed  by  its  matched  filter. 

It  was  demonstrated  that  the  Kalman  filter  can  outperform  the  RLS  algorithm, 
at  least  for  the  parameters  considered  and  the  assumputions  made.  It  was  argued  that 
the  limitation  in  performance  in  the  case  of  the  Kalman  filter  was  due  to  the  phase- noise 
and  the  residual  phase  error  caused  by  the  nonlinear  pheise  function  being  tracked.  To 
remove  the  residual  phase  error,  one  would  have  to  resort  to  a  third  order  Kalman  filter, 
thus  approaching  the  optimum  excisor  for  this  type  of  interference. 

Finally,  there  would  be  some  difference  in  performance  due  to  the  bandpcuss 
filter  at  the  receiver  input.  To  establish  how  much  would  be  related  to  the  interference 
parameters;  one  would  have  to  simulate  the  effect  of  this  matched  filter  to  obtain  precise 
results. 
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Figure  9:  Comparison  of  the  bit  error  rate  performance  between  the  RLS  algorithm  for 
A  =  0.90,  and  the  Kalman  filter  algorithm  for  d  =  0.20  rad. /sec. /volt  and  several  values 
of  low  pass  filter  bandwidth  (LPF  BW)  Bipp.  The  acronymns  NMF  and  MF  respectively 
refer  to  “No  Matched  Filter”  and  “Matched  Filter”. 
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model  designation,  trade  name,  military  project  code  name,  geographic  location  may  also  be  included.  If  possible  keywords  should  be  selected 
from  a  published  thesaurus.  e.g.  Thesaurus  of  Engineering  and  Scientific  Terms  (TEST)  and  that  thesaurus-identified.  If  it  is  not  possible  to 
select  indexing  terms  which  are  Unclassified,  the  classification  of  each  should  be  indicated  as  with  the  title.) 

LEAST  SQUARES 
RECURSIVE  LEAST  SQUARES 
KALMAN  FILTER 
SPREAD  SPECTRUM 
INTERFERENCE  SUPPRESSION 
EXCISION 
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